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In this paper we propose a scheme which allows one to find all possible exponential solu¬ 
tions of special class - non-constant volume solutions - in Lovelock gravity in arbitrary num¬ 
ber of dimensions and with arbitrate combinations of Lovelock terms. We apply this scheme 
to (6+1)- and (7+l)-dimensional flat anisotropic cosmologies in Einstein-Gauss-Bonnet and 
third-order Lovelock gravity to demonstrate how our scheme does work. In course of this 
demonstration we derive all possible solutions in (6+1) and (7+1) dimensions and compare 
solutions and their abundance between cases with different Lovelock terms present. As a 
special but more “physical” case we consider spaces which allow three-dimensional isotropic 
subspace for they could be viewed as examples of compactification schemes. Our results 
suggest that the same solution with three-dimensional isotropic subspace is more “probable” 
to occur in the model with most possible Lovelock terms taken into account, which could 
be used as kind of anthropic argument for consideration of Lovelock and other higher-order 
gravity models in multidimensional cosmologies. 
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I. INTRODUCTION 

Multidimensional paradigm being rather popular and motivated by a number of different approaches 
(from string theory to anthropic principle based on very specific and useful for life properties of gravitation 
in 3 dimensions) have also well known difficulties. We need to hide extra dimensions, so that it is reasonable 
to expect that they are contracting (or at least have been contracting during some period of Universe 
evolution), from the other hand three large dimensions are expanding almost isotropically. Such combination 
of expanding and contracting dimensions can be achieved in theories with higher order curvature corrections, 
however, the nature of this qualitative difference between large and small dimensions is often completely 
obscure. Many papers even start from decomposition of metrics as a product of “external” isotropic and 
“internal” spaces thus postulating this difference without any attempt to shed light on its origin. 

That is why any situation in which this division is not postulated a priori , but appeared due to some 


underlying principle is of particular interest. One of such situation have been found recently in Gauss- 
Bonnet gravity. Initially the solutions in question appeared when generalisation of Kasner solution for a flat 
anisotropic Universe have been studied in the regime when the Gauss-Bonnet term is dominated . It was 
found that apart from a solution in which scale factor have a power-law behavior [1, 2] (a direct analog of 
the Kasner solution) a solution with exponential time dependence of the scale factor exists [3]. It has no 
analogue in General Relativity. Later this class solution have been generalized to full Einstein-Gauss-Bonnet 
(EGB) theory (where we do not neglect the Einstein term in the action) as well as to the theory with a 
cosmological constant [4, 5]. In a more general setup such solutions (dubbed anisotropic inflation) have been 
shown to exist in theories with Ricci square corrections [6, 7] , though being absent in f(R) gravity, where 
only Gauss curvature enters into the action. 

An analysis of existence of such solutions in EGB gravity reveals an interesting fact that they exists only 
if the space has isotropic subspaces [8] (though not necessary 3-dim isotropic subspaces). The only exception 
is so-called constant-volume solutions [9], but they are special class and we do not consider them in this 
paper. Though applicability of these solutions to a realistic description of our Universe still needs more 
efforts in order to incorporate matter and exit from early time inflation, the fact that isotropic subspaces 
appears as a condition for this solution to exist is rather promising. 

One of the first attempts to find an exact static solutions with metric being a cross product of a (3+1)- 
dimensional manifold times a constant curvature “inner space”, also known as “spontaneous compactifi- 
cation”, were done in [10], but with four dimensional Lorentzian factor being actually Minkowski (the 
generalization for a constant curvature Lorentzian manifold was done in [1]). In the cosmological context it 
could be useful to consider Friedman-Robertson-Walker as a manifold for (3+1) section; this situation with 
constant sized extra dimensions was considered in [11], There it was explicitly demonstrated that to have 
more realistic model one needs to consider the dynamical evolution of the extra dimensional scale factor 
as well. In the context of exact solutions such an attempt was done in [12] where both the (3+1) and 
the extra dimensional scale factors where exponential functions. Solutions with exponentially increasing 
(3+l)-dimensional scale factor and exponentially shrinking extra dimensional scale factor were described. 

Of recent attempts to build a successful compactification particularly relevant are [13] where the dynam¬ 
ical compactification of (5+1) EGB model was considered, [14, 15], with different metric ansatz for scale 
factors corresponding to (3+1)- and extra dimensional parts, and [16, 17] where general (e.g. without any 
ansatz) scale factors and curved manifolds were considered. 

In [1] the structure of the equations of motion for Lovelock theories for various types of solutions has 
been studied. It was stressed that the Lambda term in the action is actually not a cosmological constant as 
it does not give the curvature scale of a maximally symmetric manifold. In the same paper the equations 
of motion for compactification with both time dependent scale factors were written for arbitrary Lovelock 
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order in the special case that both factors are flat. The results of [1] were reanalyzed for the special case of 
10 space-time dimensions in [18]. In [19] the existence of dynamical compactihcation solutions was studied 
with the use of Hamiltonian formalism. 

Usually when dealing with cosmological solutions in EGB or more general Lovelock gravity [20] one 
imposes a certain ansatz on the metric. Two most used (and so well-studied) are power-law and exponential 
ansatz. The former of them could be linked to Friedman (or Kasner) stage while the latter - to inflation. 
Power-law solutions were intensively studied some time ago [1, 2] and recently [3, 4, 21-23] which leads 
to almost complete their description (see also [24] for useful comments regarding physical branches of the 
solutions). Exponential solutions, on the other hand, for some reason are less studied but due to their 
“exponentiality” could compactify extra dimensions much faster and more reliably. Our first study of 
exponential solutions [5] demonstrate their potential and so we studied exponential solutions in EGB gravity 
full-scale. We described models with both variable [ 8 ] and constant [9] volume and developed general 
solution-building scheme for EGB. And now we are generalizing this scheme for general Lovelock gravity. 

In the present paper we describe a general scheme which allows us to get such solution in an arbitrary 
Lovelock gravity with arbitrary number of dimensions. This scheme is illustrated by third order Lovelock 
theory in (6+1) and (7+1) dimensions. 

As one of the objectives we want to compare solutions in the same dimensionality but with different 
Lovelock terms taken into consideration. To be specific, we will compare ( 6 +l)-dimensional solutions in 
EGB with ( 6 +l)-dimensional solutions in L± + +2 + L 3 ; then the same will be performed with (7+1)- 
dimensional solutions. The reasoning behind this comparison is simple - originally, it was EGB gravity which 
was motivated by the string theory, and third-order curvature correction which comes from string theory 
consideration, do not coincide with third Lovelock term and, through that, are not ghost-free (see [25, 26]). 
With this in mind general Lovelock theory is worse-motivated then EGB, but in higher dimensions, if we 
want to “estimate” the influence of higher-order terms, Lovelock theory could give us some “insight”. Yet, if 
we want to formally follow what comes from M/string theory, we need to stay with EGB - with both these 
arguments at hand, we will consider them both and compare the solutions one could get from both of them. 

Another important and more “physical” task is to explore the abundance of the solutions with three- 
dimensional isotropic subspace. Indeed, when dealing with higher-dimensional cosmological models, one 
needs to keep in mind that we observe only three spatial dimensions. This way we pay special attention to 
spatial splitting which have three-dimensional isotropic subspace - if these three dimensions expand while 
the remaining directions contract, this would be successful dynamical compactihcation scheme in action. 

The structure of the manuscript is as follows: first we introduce the most general equations we are dealing 
with. Then we develop the scheme for finding all possible spatial splittings in any number of dimensions 
and with any possible Lovelock terms and their combinations. After that in Section IV we apply the 


3 


scheme for (7+l)-dimensional space-times for L\ + L 2 + L 3 case to retrieve all possible spatial splittings. 
Then in Section V we describe all possible splittings in (6+1) and (7+1) dimensions and the corresponding 
solutions. After that we separate solutions which allow three-dimensional isotropic subspaces. Finally we 
draw conclusions and discuss obtained results. 


II. THE SET-UP. 


Let us spell out the conventions that we will use throughout this work. We choose to use units such that 
speed of light and gravitational constant are equal to 1; Greek indices run from 0 to D, while Latin one from 
1 to D unless otherwise stated; we also will use Einstein summation convention. 

Let us consider (D + l)-dimensional fiat space-time with Lovelock gravity. The gravitational action is 
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where n 2 is the (D + l)-dimensional gravitational constant, C m is the Lagrangian of a matter, g is the 
determinant of the metric tensor, stands for the components of the Riemann tensor, c n are constants, 
yj corresponds to the integer part of y; the generalized Kronecker delta is defined as: 

... 
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We choose a reference system in such a way that metric has the following form: 

ds 2 = - dt 2 + e 2afeW dx 2 k 
k 

Hereafter we will write a^, dk,dk instead of ak{t),dk(t),dk(t) for brevity. It is easily shown that 
R-ol = di + d 2 , R?jljl = d n d n , ji < j 2 , R[*£ = 0, {a, /3} + {/x, v}, 
the dot denotes derivative w.r.t. t. So, arbitrary component of the Riemann tensor takes the form: 

flC-fe (a. + 4) } - 

{ k i<j ) 

square brackets denote the antisymmetric part on the indicated indices. It can be shown that 
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The last term in the rhs of Eq. (6) is a total time derivative of some function; this term does not contribute 
to an equation of motion and can be omitted. Let us denote 

2 n -t _ d 


KT d = 2 n (S n - 2?r 2 S n _ 1 


E I [ V ■ £ ’” 0d = TSZ \Zi»i E c » K “ 0d P) 

h<-<32n r=l n= 0 

The superscript ’mod’ means ’modified’. Using modified Lagrangian (7) instead of initial one and varying 
the action (1) we obtain dynamical equations 


d ( 2n—2 2n 

E Cn l E + 4) E II “jr + (2n - 1) Y II “Jr \ = ^ 
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m 
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and constraint 


d 2 n 

E(2"-l)Cn E n«> = ^ 0 °, 

ji<-<j2n r= 1 


( 8 ) 
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where = c n 2 n (S n — 2n 2 G n _i). In the following we consider isotropic perfect fluid with the equation of 
state p = up as a matter source, so the energy-momentum tensor takes the form 


T 0 ° = -p, Tl = ... = TE=p 


and seek for exponential solutions such that 


ds 2 = — dt. 2 + E e 2Hkt dxf., H k = const 

k 


( 10 ) 


(ii) 


Using the notations of the Sec. II we see that ak(t) = Hkt, = const. Substituting it in (8)-(9) we see 
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It could be useful to rewrite Eqs. (12)-(13) in terms of elementary symmetric polynomials (see e.g. [27]): 

(14) 

In what follows we will also use helpful notation for e n with parameters H / C1 ,..., excluded: 

72 
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It is easy to check that 
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Then with (14)—(16) taken into account Eqs. (12)—(13) read as 

d d 

0 = £ i = Cn ( e l e 2n-l ~ e 2n) “ WK, i = 1, . . . , D ] ( 2n “ l)Cne 2 n = (17) 

n=l n =1 

Since we assume Hi = const, it follows from Eqs. (12)-(13) that p = const, so that the continuity equation 

P+{P + P)J2 H ^ = 0 ( 18 ) 

i 

reduces to 


{p + p)Yl Hi = 0 ’ ( 19 ) 

i 

which allows several different cases: a) p = 0 (vacuum case), b) p + p = 0 (A-term case), c) Yi^i = 0 
(constant volume case which we call CVS for brevity) and their combinations: d) vacuum CVS and e) 
A-term CVS. So, all possible exponential solutions can be divided into two large groups: solutions with 
constant volume and solutions with volume changing in time (non-constant volume solutions); for the latter 
case we have only two possibilities: vacuum case and A-term case; on the contrary, the first case does not 
impose constraints on choice of matter a-priory. 

Let JV£N and let [ N ] be the set of the first N natural numbers: 

[N] = {l e N| l < N} (20) 


Initial system of D dynamical equations is equivalent to a system that consist of D — 1 difference equations 
and one dynamical equation: 

V i € [D] (£i = 0^ V i € [D — 1] (£i+ 1 — £i = 0^ A 3m € [D] ( '£ m = 0^ (21) 


It is easy to check that 
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Let us introduce the following notations: 
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Then we have 


( 22 ) 


(23) 


£i+i £i — ei Cj j +1 Ci.i+i — 0 <—> C i i+1 — 0 V — 0 V ^ Hk — 0 (24) 

k 

Equation Yk = 0 leads to so called constant volume solution (see [8, 9] for details) and requires separate 
consideration, so in what follows we deal with equations C\ i+ , = 0 and Cf i+1 = 0 only. 
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One more thing needs to be explained before continuing with the general scheme, namely, the coupling 
constants ( n . In previous papers dedicated to study of the exact solutions in EGB and Lovelock gravity [4, 
5, 8 , 9, 22, 24] we used different couplings constant and for a reader’s convenience in comparing solutions 
we write down their relation (remember, ( n = c n 2 n (S n — 2n 2 G n _i) with S n given in ( 6 )): 

Ci = -2,(2 = -8a, C 3 = -144/3. (25) 

In the absence of (2 and (3 natural normalization for dimensional units would be 

- (i = 8itG => 4t tG = 1, (26) 

so that when we put some dimensional numerical values, we do it with (26) in mind. 


III. THE GENERAL SCHEME. 


In this section we describe the general scheme that allows to obtain a wide class of solutions with isotropic 
subspaces in the framework of general Lovelock model in (D + 1) dimensions. 

Let m G Z>o, N G N; m < IV. Let us introduce the following definitions: 


1. Let ZpL be a collection of all m-element subsets of the set [TV]: 

l$ ] = {zc[N]\\Z\ = m} 

Notation \Z\ in (27) stands for cardinality of a set Z. 

2. Let r, k r . l r _ 1 € N, Iq = U, X r G Xjx _ 1 ; let us introduce the following notations: 

Jo = [D ], Jk = Jk-i\Jk-i, Jk-i = Jk-\ u Jk-u where 


(27) 


Jk -1 = {j 
Jk -1 = {is+A 


js € Jk—i A sG A), A s < —^— j | 
js+k € 77fc_i A s G ^ A s > -^5 — J | 


(28) 


In what follows we use special notation y for elements of i.e. if we write y J somewhere it means 
there exists a set Jk such that y% G 77fc. 

1st step. For a given solution some of equations £i+i — £\ = 0 are satisfied due to C / i+1 = 0, others 
are satisfied due to C 2 i+1 = 0, thus the initial system becomes equivalent to a system that consist of several 
(say, k\) equations C} i+l = 0, several (say, l\) equations Cf i+1 = 0 and one dynamical equation: 

3 k\ G Z>o 3 X\ G Zp_ 1 j V i G A) (C/,i+i = 0^ 

3lrG Z> 0 3 y x G Vj G ^ (c 2 j+1 = o) 

AinLi = 0 a x 1 \jy 1 = [D-i\ ( 29 ) 

£\ = 0 

+ l\ = D — 1 


Mi G [D] ($ = 0) 
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Note that if k\ = 0 (Zi = 0) then X\ = 0 (34 = 0 ). 

2nd step. The idea is to consider equations C 2 J+1 = 0, j G 34 as a new basic equations, hnd the 
respective difference equations and obtain result analogous to (29). The difference C\ ■ — C 2 ■ is factorized 
iff one of the elements of the pair (j' 3 , j' 4 ) equals to one of the elements of the pair (ji, J 2 ); for example, let 
us assume that j '4 = j 2 , then 


r>2 _ n2 _ nl n2 

1J2 _ U jl ,j3 U il,j2j3 


(30) 


According to (23) C\, i/+l = 0 H = Hjj + \, therefore we can identify indices i',i' + 1; it is easy to 

check that one can always find such f 6 T, j' G 34 that %' = j' + 1 (or «' + 1 = j')-, using these facts we 
replace the set 34 in (29) by J\ . By analogy with (21) we have 
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A 2 n 34 — 0 a X2 u 34 — [4 — 1 ] 
fe 2 + / 2 = 4 — 1 


(33) 


Taking into account 1st and 2nd steps we see that 

Cli+l = 0 for all i G Xi G 

C ^ +2 =°> fora llPG^GI ^j 

C ^ +1 ^ +2 =°’ forall^G^Gl^ 

C 2 1 1 =0 for some A:q G [Zi] 

fc o fc o +1 

£ m = 0 for some m G [D\ 
k\ + k~2 + Z 2 = D — 2 


f£i = 0 
[ £d = 0 


(34) 
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r-th step. Continuing this procedure at r-th step we have 



C pi,pi+i = 0 for all -Pi € x i € X [b- 1 ] 
C\ i =0 for all p 2 £ X 2 £ , 

V P2’ V P2 +2 11 J 


C r _j. = 0 for all p r € X r € Z(? r n 

C 2 r _i r _i = 0, for all q r E TV E , _j] ■ lr > 1 


Vq r ? ■ ■ ■ )U q r -|~7” 

C 2 i i = 0 for some /cq € [Zi] 


(35) 


r 2 

^ r—1 r—1 

y r— 1 ’*••’ 2 /, r —1 , T 

fc o h> + 7 - 1 

£ m = 0 for some m E [D] 

ki “h ... “I - k r + l r = D — v 


= 0 for some 1 E [Z r _i] 


The process terminates when l r < 1 for some r = r* E N; the initial system becomes equivalent to the 
system consisting of k\ +... + k r * = D — r* equalities of Hubble parameters (equations C 1 i = 0) and r* 

y Pi ’ y Pi +2 

additional conditions. The maximal value of r* is (2d — 1). Indeed, when r* = r^ ax = (2d — 1) degree of the 
last symmetric polynomial in the sum (23) becomes zero and C 2 • t = Q- Equation C 2 = Q = 0 


means vanishing of the highest Lovelock term in the Lagrangian and must be omitted, in other words one 
must set l r ^ ax = 0, k r ^ ax = l r ^ ax -1 — 1- Thus, r* can vary from 1 to r^ ax and, respectively, (&q + ... + k r *) 
varies from (D — 1) to 1 for even D and to 2 for odd D. 

Each collection of (k\ + ... + k r *) equalities of Hubble parameters gives a number of splittings of space 
into isotropic subspaces. Splittings corresponding to various values of (k\ + ... + k r *) are represented in 
the Table I. We see that even-dimensional spaces has one more splitting as compared with odd-dimensional 
one. Numbers in the column "Splitting" means numbers of equal Hubble parameters; braces stand for 
pairing of Hubble parameters that give rise to the next splittings; subscripts after round brackets are used 
to indicate the number of units in these brackets. For example, record 3 + 2 + (1 +1 + ... + 1 )d -5 means 
that {iZi,..., Ho} = {H, H, H,h,h, Hq, ..., Hp>}, H§ + ... + Hp> and at the next step one can obtain the 
following splittings: 


5 + (1 + ... + 1)_D—5, 3 + 3 + (1 + ... + l)z?—6, 3 + 2 + 2 +(1 + ... + 1 )d_7 


(36) 
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k\ + ... + k r * 

Splitting 

Number of additional conditions 

1 

2 +jT+T+... + l) 2,_2 

D-l 

2 

3 + (1 + 1 + . . . + 1)d—3 

2 + 2 (1 + lyf . . . + 1 ) D —4 

D-2 

3 

4 + (1 + 1 + . . . + 1) D —4 

3 + 2 (1 . . . + 1) D —5 

2 ^ 2+2 + (1 -|- 1 -|~ ... + 1 )d —6 

D-3 




D-l 

isotropic 

1 


Table I. Several splittings of D-dimensional space. 


IV. APPLICATION OF THE GENERAL SCHEME. 


In this section we consider application of the of the scheme described above to the theory with three 


Lovelock terms in (7 + 1) dimensions. Let 


Vh2 = 


Cl 

c 2 ’ 



(37) 


I. r* = r^ ax = 5; Z 5 = 0, k\ + ... + k§ = 2. It is easy to check that 


k\ + ... + &5 — 2 

t 

1. (k 5 = 2,/ 5 = 0); (k A = (1. l 4 = 3); ( k 3 = 0 ,l 3 = 4); (k 2 = 0,l 2 = 5); (h = 0, h = 6), or 

2. (fc 5 = 1, h = 0); (h = 0 ,U = 2); (k 3 = 0,l 3 = 3); (k 2 = 0 ,l 2 = 4); {h = 1, h = 5), or 

3. (k 5 = 1 ,k = 0); (k 4 = 0 ,l 4 = 2); ( k 3 = 0J 3 = 3); (k 2 = 1 ,l 2 = 4); {h = 0, h = 6), or 

4. (k 5 = l,l 5 = 0); (k 4 = (1. l 4 = 2); ( k 3 = 1 ,l 3 = 3); (k 2 = 0, l 2 = 5); (fci = 0, h = 6), or 

5. (h = l,h = 0); (M = U4 = 2); (k 3 = 0,l 3 = 4); (k 2 = 0 ,l 2 = 5); (h = 0 ,h = 6) 


(38) 


We consider subcases (1) and (2) in more details; other subcases can be considered analogously. 


(1) It is easy to check that 


X! = 0 , 34 


[ 6 ], *2 


0 , 34 = [5], X 3 = 0 , 34 
Ji = J2 = Jz 


[4] , *4 = 0, 34 = [3] , X 5 = {1,2} , 34 = 0 
Ja = [7] 
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fel + . . . + k r * 

Splitting 

Number of additional conditions 

2 

3+1+1+1+1 

2+2+1+1+1 

5 

3 

(4+1 + 1 + 1) 

(3 + 2 +1 + 1) 

(2 + 2+ 2 + 1) 

4 

4 

(5 + 1 + 1), (3 + 2 + 2) 

(3 + 3 + 1), (4 + 2 + 1) 

3 

5 

(6 + 1), (5 + 2), (4 + 3) 

2 

6 

isotropic 

1 


Table II. All possible splittings of 7D space. 


Without loss of generality we can choose fog = ... = = m = 1, then 


f £1 = 0 
1 £7 = 0 


£•16 = £"27 = 0 
£"12345 = 0 
£"1234 = 0 
£"123 = 0 

£•12 = 0 

£l=0 


' H x = H fi = x, H 2 = H 7 = y 
H 3 = z, H 4 = u, H 5 = v 
x = -y, C2 - C3Z 2 = 0 
C 2 CU + v) - ( 3 (u + v)x 2 = 0 

Ci + (2(zu + zv + uv — x 2 ) — C3 (zu + zv + uv)x 2 = 0 
3 

E (n (+L-1 - e\ n ) = ujx (*) 

n= 1 


(40) 


We see that (40) gives us (2+ 2 + 1 + 1 +1) splitting and necessary conditions for a set (x, x, y, y, z, u, v ) to 


be a solution of system under consideration: 


x 2 = + ’ ^32 > 0, 1^12^32 = 1, y = -x; z, u, v G R (41) 

It is easy to check that the same relations are obtained for any other choice of fcg, ..., Aiq; to get final solution 

one should substitute these relations into dynamical equation (*) and constraint (13). (2) We have 

X x € Xl let X\ = {1} , then 33 = {2,..., 6} , Jx = {2,..., 7} 

*2 = 0, 33 = [4] , X 3 = 0 , y 3 = [3] , x 4 = 0 , y 4 = [ 2 ], *5 = {1}, 3% = 0 (42) 

= Jz = J 2 = Ji = {2, ■ ■ •, 7} 

Without loss of generality we can choose fcg = ... = = m = 1, then 

(£1 = 0 r c} 2 =c \ 7 = 0 r ^ = h 2 = h 7 = x 

< ... 4—b | C23 = £"234 = £-2345 = £"23456 = 0 < ' == ’ > j -^3 = Vi H 4 = Z , H§ = U, H§ = Z, (43) 

{ £7 = 0 [ Ex = 0 1 X = 0, C2 = 0 

Equality C2 = 0 implies that the second Lovelock term is absent which contradicts to formulation of the 
problem (see Sec. II). So, in the case under consideration splitting (3 + 1 + 1 + 1 + 1) is not implemented. 
It is not difficult to verify that the same result is obtained for any other choice of X\ and /eg,..., fcg. One 
can easily make sure that the rest subcases (3)-(5) does not give us any new results. 


11 













II. r* = 4 . In this case k\ +.. . + k 4 + l 4 = 3; there are a lot of combinations like (38) which give such sum, 
so we consider only three of them; one can easily check that another combinations give the same results. 
(1) (k\ = 0, l\ = 6), (k /2 = 0, h = 5), (k 3 = 0, £3 = 4), (k 4 = 3, U = 0). It easy to check that 


xi = 0, y 1 = [e], *2 = 0, y 2 = [5], 

Ji = j 2 - 

Without loss of generality we can choose k^ = ... - 


X 3 = 0, 0>3 = [4] , 
Js = Ja = [7] 

fcjj = m = 1, then 


*4 = {1,2,3}, 314 = 0 


(44) 


f Si =0 

l ^7 = 0 


cj 5 = 0 

^26 = £37 = 0 

^1234 = 0 
£"123 = 0 

Cj 2 = 0 

£l=0 


H\ = H Fj = H 2 = Hq = y, H 3 = H 7 = z, H 4 = u 
C 2 + C 3 (xy + xz + yz) = 0 

( 2 (% + y + z + u) + (3 (u[xy + xz + yz\ + xyz ) = 0 

Cl + C 2 (u[x + y + 2 z\ + xy + 2 xz + 2 yz + z 2 )+ 

+( 3 (u[ 2 xyz + xz 2 + yz 2 ] + xyz 2 ) = 0 

3 

E Cn (eje^n-! - e\n) = UK (**) 

n= 1 


(45) 


We see that (45) gives us (2 + 2 + 2 + 1) splitting and necessary conditions for a set (x, x, y, y, z, z, u) to be 
a solution of system under consideration: 


X 2 = ^ 32 \ V>32 > 0, "012^32 = 1, y = -X, Z,U<ER 


(46) 


It is not difficult to verify that the same result is obtained for any other choice of kg,kg; to get final 
solution one should substitute relations (46) into dynamical equation (**) and constraint (13). 

(2) (ki = 0, l\ = 6), {k -2 = 1, h = 4), (k 3 = 0, l 3 = 3), (k 4 = 2, Z 4 = 0). We have 


X\ = 0, 33 = [6], Ji = [7] 

*2 GlJ, let A7 2 = {1} , then 3^2 = {2,..., 5} , J 2 = {2,...,7} (47) 

*3 = 0, 3^3 = [3], J 3 = J 2 - X 4 {1, 2} , 3^4 = 0 


Without loss of generality we can choose ky = ... = = m = 1, then 


f £1 = 0 
l ^7 = 0 


C43 = 0 

^26 = £37 = 0 

^2345 = 0 
^234 = 0 

C 2 3 = 0 

£l=0 


Hi = H 3 = H 7 = X, H 2 = H 6 = y, H 4 = z, H 5 = u 
C 2 + (3(x 2 + 2 xy) = 0 

C 2 ( 2 x + y + u) + ( 3 (n[x 2 + 2 xy] + x 2 y) = 0 

Ci + C 2 ([2x + y](z + u) + x 2 + 2 xy + zu) + 

+C 3 ( zu[x 2 + 2 xy] + (z + u)x 2 y) = 0 
3 

E Cn (eje^ - el n ) = ux (* * *) 

n= 1 


(48) 


We see that (48) gives us (3 + 2 + 1 + 1) splitting and necessary conditions for a set ( x , x, x, y, y, z, u ) to be 
a solution of system under consideration: 


V 2 = ^32 > 0, ^ 12^32 = 1, x = -y, z,u£ R 


(49) 
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It is not difficult to verify that the same result is obtained for any other choice of X 2 and k $to get 
final solution one should substitute relations (49) into dynamical equation (* * *) and constraint (13). 

(3) ( k\ = 2, li = 4), (k -2 = 0, l 2 = 3), (k 3 = 0 ,l 3 = 2), (/c 4 = 1, l 4 = 0). We have 


*1 € X|, let * 4 = {1, 2} , then 33 = {3,4,5,6} , J x = {3,..., 7} 

X 2 = 0, 3^2 = [3], J 2 = J 1 , ^3 = 0, 3^3 = [2], Jz = J 2 = Jl, X 4 {1} , 3^4 = 0 


(50) 


Without loss of generality we can choose k^ = ... = k$ = m = 1, then 


f £1 = 0 

[ ^7=0 


£"12 = ^23 = ^37 = 0 
^3456 = 0 

£345 = 0 
C 3 2 4 = 0 
£\ = 0 


' Hl = H 2 = H 3 = Hj = x , 

H 4 = y,H 5 = z,H 6 = u 
C 2 + 3C3^ 2 = 0 

< C2(3x + u) + ( 3 ( 3 ux 2 + x 3 ) = 0 

Ci + C 2 (3x 2 + 3 x(z + u) + zu) + 

+C 3 (3 x 2 zu + x 3 (z + u )) = 0 
3 

E Cn ( e l e 2n-i - e L) = UX 


(51) 


Equality C2 = 0 implies that the second Lovelock term is absent which contradicts to formulation of the 
problem (see Sec. II). So, in the case under consideration splitting (4 + 1 + 1 + 1) is not implemented. It is 
not difficult to verify that the same result is obtained for any other choice of X\ and k^..... . 

Cases r* = 1,2,3 can be considered analogously. We do not give detailed computations due to their 
awkwardness, but in the next section we describe all possible splittings in detail for (7 + 1) and (6 + 1) 
dimensions. 


V. NON-CONSTANT VOLUME SOLUTIONS IN (6 + 1) AND (7 + 1) DIMENSIONS. 

Below we give non-constant volume solutions in an implicit or an explicit (wherever it is possible) form 
for the models with three Lovelock’s term (L\ + L 2 + L3) and two Lovelock’s term (Li + L 2 ) in (7 + 1) and 
(6 + 1) dimensions. 

We have two free parameters in EGB model: -012 and x. Let O2 be a two dimensional space of parameters 
(ijj 12,«); let 0| o1 be a subspace of space O2 such that solutions exist iff ('012,x) £ f2| o1 , dim (^1°') ^ 2; we 
point out explicitly the dimensionality of subspace n| o1 if dim < 2 and highlight the corresponding 

splittings by boxes. 

In the model with three Lovelock’s term there exist three free parameters: vpi 2 ,x,' l p 32 - In a similar way 
we denote via ^3 three-dimensional space of parameters (ipi 2 , ^32) and introduce S3| o1 - subspace of space 
H3 such that solutions exist iff (0i2, ^32) € fl| o1 , dim < 3; as in the case of EGB model we point out 

explicitly the dimensionality of subspace n| o1 if dim < 3 and highlight the corresponding splittings 

by boxes. 
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A. (7 + l)-dimensional spacetime. 


1. (2 — 2 + 1 + 1 + 1) = (x, x, —x — x, y, z, u ) 


L\ + L2 + L3 

necessary conditions: 

X 2 = ^32 > ^32 > 0, 0 12032 = 1, 2/, U € R 

substitution into equation of motion and constraint leads to 
w = -1, * = ''A-7 2 \ dim = 1 

This case involve two subcases: 

• if y = z then (2 — 2+1 + 1 + 1) is transformed into 

(2 — 2 + 2 + 1) = (x, x, — x, —x, z, z, u ) 

• if y = x then (2 — 2 + 1 + 1 +1) is transformed into 

(3 — 2 + 1 + 1) = (x, x, x, —x, —x, z, u ) 


L\ + L-2 


No solutions 


(52) 


2. (3 — 3 + 1) = (x , x, x, —x , —x, —x, y) 


1 


Let us denote A± =- 1 ± \J\ — 012032 

032 1 


(53) 


x 2 = 


L\ + L 2 + L 3 


necessary conditions: 
A + , 012032 <1, 032 > 0 
A-, 012032 < 1, 012 > 0 


2/e 


Substitution into equation of motion and constraint 
leads to 

C 3 X 6 - 3Cix 2 


a) oj = — 1, « = 


2 n dim(ll| ol )=2 
b) x = 0, 0 i 2 032 = dim = 1 


L] + L2 


necessary conditions: 

u = -1, x 2 = -yC '012 > 0 , y G 

Substitution into equation of 
motion and constraint leads to 

3 Ci 2 


= — 


2C2 


Ci, C2 < 0, dirndl 01 ) 


= 1 


(54) 
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3. (3 + 3 + 1) = (x, x, x, y, y, y, z) 


L i + L 2 + 

necessary conditions: 

, , 1~ 012032 ,-1 

(x + yY =- — -, xy = —032 , ^ € R 

(012032 < 1 A 032 > 0) V (012032 > 5 A 0 32 < 0) 


Substitution into equation of motion and constraint leads to 
Ci(3x 4 + 9 x 3 y + llx 2 y 2 + 9 xy 3 + 3 y 4 ) 


a) oj = — 1, x = — - 


b) x = 0, let £ = — then 

x 


x 2 + 3 xy + y 2 


dim (n| o1 ) 


= 2 


3£ 4 + 9£ 3 + 11£ 2 + 9£ + 3 = 0, & « -1.662, £ 2 ~ -0.602 

V = dim ffl! 01 ) = 1 


x 2 = — - 


£032 ' 


L\ + L-2 


No solutions 


(55) 


4. (3 + 2 + 2) = (x, x, x, y, y, z, z) 


Let us denote 

B\{x) = 0f 2 X 8 + (012032 - + (7032 - 012032) a;4 + (012032 ~ 3012032 “ 2)x 2 + 0 12 

.B 2 (x,y) = — Ci(9x 4 + 14x 3 y + llx 2 y 2 + 4xy 3 + y 4 ) 

5 3 (£) = 1026£ 5 + 6012£ 4 + 15078£ 3 + 22159£ 2 + 16524£ + 5196 ( 56 ) 

1?4(x, y, z) = 3x 5 (y + z) + 6x 4 (y + z) 2 + 4x 3 (y 3 + 4 y 2 z + 4yz 2 + z 3 )+ 

2x 2 (y 4 + 5 y 3 z + 8 y 2 z 2 + 5yz 3 + z 4 ) + 4x(y 4 z + 2y 3 z 2 + 2y 2 z 3 + yz 4 ) + 2y 4 z 2 + 2y 3 z 3 + 2y 2 z 4 

B$(x, y, z) = x 3 y + x 3 z + 2x 2 y 2 + 4x 2 yz + 2x 2 z 2 
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L\ + L 2 + L3 

necessary conditions: 


~ _ (012032 ~1)®± \JB\{x) _ 0 32 X 2 z + 2 x + z 

( 0 32 ® 2 - l) 2 ’ V 032®(® + 2z) + 1 

with x such that -Bi(x) ^ 0 

Substitution into equation of motion and constraint 
leads to 
2 B 4 (x,y,z) 


a) oj = — 1, x = 


B 5 (x,y,z) 


b) x = 0, let £ = —, x = — then equation B 4 (x,y,z) = 0 

x ' x 

take the following form: 

2 (£ + l)V + 2 (£ + 2) (£ + l)V + 2 ( e 2 + 2 £ + 3) (£ + i) 2 x 2 
+ (4£ 4 + 10£ 3 + 16£ 2 + 12£ + 3)x + 2£ 4 + 4£ 3 + 6£ 2 + 3£ = 0 


Special solution: 


012 


x =j^b 3 (Z), y = £®, 2 = -' 


2£x 


432 „ * , 3 ^ + 2 
18£ 6 + 120£ 5 + 350£ 4 + 603£ 3 + 606£ 2 + 324£ + 72 = 0, 
6 « -2.556, £2 » -0.902 


Li + L'2 


necessary conditions: 


z = —x± \Jip 12 - 2x 2 , 
y = —2x — z 

x such that x 2 < ip 12 > 0 

Substitution into equation of 
motion and constraint leads to 


W = —1, H = 


B- 2 (x,y) 


3x 2 + 2xy + y 2 


5. (4 + 2 + 1) = (x, x, x, x, y, y, z) 


L\ + L-2 + L 3 


necessary conditions: 

032®® - 3 lp 32 X 4 - 3 ('012'032 - 2) x 2 - 012 = 0 
x(0 32 x 2 + 3) 


y = 


z G 


3032® 2 + 1 

Substitution into equation of motion and constraint 
leads to 

Ci(15x 4 + 45x 3 y + 48x 2 y 2 + 24xr/ 3 + 8y 4 ) 


w = —1, x = — - 


3x 2 + 9xy + 8 y 2 
dim = 2 


L\ + L 2 


necessary conditions: 

® 2 = 012 > 0 
6 

y = —3x, z€l 

Substitution into equation of 
motion and constraint leads to 

13Ci 2 


LO = —1, X = — 


6(2 

dim (fif= 1 


Cl, C2 < 0 


(58) 
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6. (4 + 3) = (x,x,x,x,y,y,y) 


Let us denote 

F\(x) = 2 X 6 - 3^32X 4 - 3(^12^32 - 2)x 2 - Ip 12 

F 2 (x,y ) = — C2(6 x 5 + 36x 4 y + 66 x 3 y 2 + 51x 2 y 3 + 16xy 4 ) — £i(12x 3 + 27x 2 y + 22 xy 2 + 9y 3 ) 
F 3 (x,y) = —3Ci(5x 4 + 12x 3 y + llx 2 y 2 + 6xy 3 + y 4 ) 


(59) 


L\ + L 2 + L3 


L\ + L2 


necessary conditions: 

-V> 3 2X 3 - 3x ± , , 1 , „ / w n 

y =-x-, x such that L+x) ^ U 

3+32 x^ + 1 

Substitution into equation of motion 
and constraint leads to 

X , F 2 {x,y) 

a) 2^T3i/ 

b) « = 0, F 2 {x,y) = 0 


necessary conditions: 

3x 2 + 6 xy + y 2 = —ip\ 2 

Substitution into equation of motion 
and constraint leads to 

a) u = -l, x =-ip^F 3 (x,y) 

b) x = 0, let £ = — then 

x 

{ £ 4 + 6£ 3 + 11C 2 + 12£ + 5 = 0, 



-3.874, £ 2 ~ -0.743 
^12 _ t 

e 2 + ee + 3’ y ^ 


(60) 


7. (5 + 1 + 1) = (x, x, x, x, x, y, z) 


L\ + L 2 + L 3 

necessary conditions: 

x 2 = ~ips 2 , tp32 < 0, ip 12^32 = 5, |/, + l 

Substitution into equation of motion and constraint leads to 
u = —1, x = — 7V ; 32 1 ) dim (R| ol j = 1 


Li + +2 


No solutions 


(61) 


8. (5 + 2) = (x,x,x,x,x,y,y) 


Let us denote 

G\(x , y) = —C2(15 x 5 + 60x 4 y + 80x 3 y 2 + 20x 2 y 3 ) — £i(10x 3 + 40x 2 y + 16xy 2 + 4y 3 ) (62) 

G 2 (x, y) = —Ci(45x 3 + 40x 2 y + 16xy 2 + 4y 3 ) 
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L\ + L2 + L3 


y = 


necessary conditions: 
032® 4 + 6x 2 + '012 


x € 


4a; (1 + 032 -t 2 ) 

Substitution into equation of motion and 
constraint leads to 

Gi(x,y) 

a) w = —1, h =-—^ 

x + Ay 

b) x = 0, Gi(x,y) = 0 


L\ + L-2 

necessary conditions: 
6x 2 + 012 


y = 


4x 


X £ 


Substitution into equation of motion and 
constraint leads to 
G 2 (x,y) 


a) co = — 1 , x = 


2 (3.x + 2 y) 


b) x = (1. let £ = — then 

x 

4£ 3 + 16£ 2 + 40£ + 45 = 0, £ « -1.870 

2 ’012 t 

2(2£ + 3) 


( 63 ) 


9. (6 + 1) = ( a; ^ a; ^ a; ^ a; ^ x ^ x ^ y ) 


Let us denote J± 


1 

032 


-1 ± 


1 _ 012032 
5 


(64) 


L 1 + L2 + H3 


x 2 = 


necessary conditions: 
J+i 012032 < 5, 012 < 0 
J_, 012032 < 5, 032 < 0 


y e 


Substitution into equation of motion and 
constraint leads to 

35(0x 6 — 210x 2 ( ol 

a) cu = — 1, ^ = --- . dim ^S2 3 

I K / x 

b) jc = 0, 012032 = -p dim ^| ol J = 1 


= 2 


Li + L2 


necessary conditions: 
x 2 = pp, 012 > 0, y € R 

Substitution into equation of motion and 
constraint leads to 

21Ci 2 


(65) 


uj = —1, x = — - 


2OC2 

dim = 1 


i C15C2 < 0 


10. (7 + 0) — ^X ^ X j Xy Xy Xy Xy X ^ 


Let us denote 


K± = 


2032 


- 1 + 01 - 


4012032 

15 


L± = 


012 

10 


- 1 ± 01 - 


20 x 


21 0120 


( 66 ) 
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L\ + L 2 + L% 

Substitution Hi = x into equation of motion 
and constraint leads to 

a) ui = —1, x = —35C3X 6 — 105C2® 4 ~ 2 \Q\x 2 

15 

K + , ^12^32 < -r, 1 P 12 < 0 

b) h = 0, x 2 = 

^12 V ; 32 < —, ^32 < 0 


Tl + -^2 


Substitution Hi = x into equation of motion 
and constraint leads to 


a) u = — 1, x 2 = 


L 


+) 




jz 21 

^12(1 ^ 20’ ^ < ° (67) 

5/ 21 

< 7^7: i i >\2 < 0 


V , i2Ci 20 


b) x = 0, x 2 = —ip 12 < 0 
5 


One can see that in (7+1) dimensions solutions for Li + +2 + T3 are more abundant then solutions for 
Li + L 2 (EGB) - probably due to the higher number of degrees of freedom - we will discuss this aspect in 
detail in the appropriate section. For the reader’s convenience we summarize all the solutions in Table III. 


Table III. Summary of (7+l)-dimensional solutions in Li + L 2 + L 3 and EGB. 


Splitting 

nonvacuum Li + L 2 + T3 

vacuum L\ + L2 + L3 

nonvacuum EGB 

vacuum EGB 

(2 — 2 + 1 + 1 +1) (52) 

yes 

no 

(3-3 + 1) (54) 

yes 

no 

(3 + 3 + 1) (55) 

yes 

no 

(3 + 2 + 2) (57) 

yes 

no 

(4 + 2 + 1) (58) 

yes 

no 

yes 

no 

(4 + 3) (60) 

yes 

(5 + 1 + 1) (61) 

yes 

no 

(5 + 2) (63) 

yes 

(6 + 1) (65) 

yes 

no 

(7 + 0) (67) 

yes 


B. (6+l)-dimensional spacetime. 


1. (2 - 2 + 1 + 1) = (x, x, —x, —x, y, z) 

L\ + +2 + T3 
necessary conditions: 
x 2 = ^32 • y,z eR, ip 12^32 = 1 
substitution into equation of motion and constraint leads to 
w = - 1, x = ip^, dim(ll| ol ) = l 


L\ + L 2 


No solutions 


( 68 ) 
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2. (2 + 2 + 2) = (x,x,y,y,z,z) 


Let us denote 


Ml(y) = 4^322/ 4 + [^12^32 - 6^12^32 - 3 \y 2 + 4^12 
M 2 (x, y) = x 3 y + x 3 z + x 2 y 2 + 2 x 2 yz + x 2 z 2 + xy 3 + 2xy 2 z + 2xyz 2 + xz 3 + y 3 z + y 2 z 2 + yz 3 


Li + L 2 + L3 

necessary conditions: 

_ (1 - ^12^32)2/± VMi(y) __ _ x + y 

2(V’32?/ 2 -l) 1 + ^32Xy 

with y such that Mi(y) > 0 

Substitution into equation of motion and constraint 


leads to 


oj = —1, % = 


M 2 (x,y) 
xy + xz + yz 


L\ + +2 


No solutions 


3. (3 + 2 + 1) = (x, x, x, y, y, z) 


Let us denote 


Ni = 3 - 12^32 ± \/(3 - V’i2'i/’32) 2 - 4V>12'032 , N 2 (x, y) = -Ci(3x 3 + 6x 2 y + 4xy 2 + 2y 3 ) 


Li + L 2 + L3 

necessary conditions: 


L\ + L 2 


7V[ + , for 


N x , for 


V>32 > 0 


V'12'032 < 1 V’12V’32 > 9 


1 + ^32X 2 

•032 < 0 


■032 > 0 , 1pi2 > 0 

^ 12^32 < 1 


^32 < 0 


^12^32 > 9 or lpi 2 > 0 


Substitution into equation of motion and constraint 


leads to 


a) a; = —1, x = 


N 2 (x,y) 
x + 2 y 


b) x = 0, let £ = — then 

x 


2£ 3 + 4£ 2 + 6£ + 3 = 0, £ « -0.722 
2 £ + 2 t 

® = - 7 - 7 —, y = £x 

£+32 


necessary conditions: 

x 2 = -p ^12 > 0 , 
y = —2x, 

Substitution into equation 
of motion and constraint (72) 
leads to 


1 ^>1 

w = -l, x = —7^, 

42 

Cl) C2 < 0, dim (Q| o1 ) = 1 
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4. (4 + 1 + 1) = (x, x, x, x, y, z 


L\ + L/2 + L 3 
necessary conditions: 

X 2 = —3lf>22 1 ^32 < 0 , ^12^32 = 9, y, z € R 

Substitution into equation of motion and constraint 
leads to 

a j = —1, ^ = — 15-032 1 , dim = 1 


L\ + L 2 


No solutions 


5. (3 + 3) = (x,x,x,y,y,y) 


Let us denote 

Pi(x) = -ip 32 x 4 + (3 - *l>i2il>32)x 2 - Vh2 
P 2 (x, y) = -C2xy(4x 2 + 7xy + 4y 2 ) - Ci(3x 2 + 4 xy + 3y 2 ) 
P 3 (x , y) = —3Ci(x 4 + 4 x 3 y + 5x 2 y 2 + 4xy 3 + y 4 ) 


Li + L 2 + L 3 

necessary conditions: 
-2X db >/Ptfx) _ . . 


y =---, x such that P±(x) ^ 0 

W32X Z + 1 

Substitution into equation of motion 
and constraint leads to 

a) u = -1, x = P 2 (x,y ) 

b) x = 0, P 2 (x, y) = 0 


L\ + L 2 

necessary conditions: 

y = —2x ± sj’&x 1 - ^12 
with x such that x 2 ^ 

O 

Substitution into equation of motion 
and constraint leads to 


a) lx = —1, = 


P3(x,y) 

x 2 + 4 xy + y 2 


b) x = 0, let f = — then 

x 


£ 4 + 4£ 3 + 5£ 2 + 4 £ + 1 = 0 , £ = -- ± 

2 ^12 t 

1 ~ - g 2 + 4£ +1 1 y = (x 
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L\ + L2 + L3 


y = 


necessary conditions: 

3x 2 + 0 i2 


x G 


03220 + 3x ’ 

Substitution into equation of motion and constraint 
leads to 

a) u = -1, x = Qi(x,y) 

b) x = 0, Qi(x,y) = 0 


L\ + L '2 

necessary conditions: 

3x 2 + 0i 2 


y = 


3x 


x G 


Substitution into equation of motion and 
constraint leads to 
Q2(x,y) 


a) u = —1, x = 


x + y 

b) x = (1. let £ = — then 

x 

£ 3 + 3£ 2 + 6£ + 5 = 0, f = -1.322 


x = — 


012 


3(C + 1) 


, y = & 


( 77 ) 


7. (5 + 1) = (x, x, x, x, x, y) 


Let us denote U± 


1 

'032 


3 ± \/9 — 012032 


(78) 


L\ + L2 + L3 
necessary conditions: 


L\ + L2 


j/€R, .t 2 


for 


032 > 0 

or 

■012 < 0 


■032 < 0 
0 < 012032 < 9 


, for 


"032 < 0 
012032 < 9 


Substitution into equation of motion 
and constraint leads to 

a) u = — 1, « = — 1502a: 4 — lOCix 2 , dim ^Q| o1 j = 2 

b) ^ = 0, 012032 = dim (^3 Cl ) = 1 


necessary conditions: 

x 2 = 012 <0, y G R 

b 

Substitution into equation of motion 
and constraint leads to 

u = — 1 , x = — ^ 12 , dim = 1 


(79) 


8 . (6 + 0 ) — (x 1 1 3:. 3:. 3:. 3:. x ) 


Let us denote V± 


1 

2032 


— 9 ± 081 — 12012032 , W± 



(80) 
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L\ + L 2 + L% 


Substitution Hi = x into equation of motion and 
constraint leads to 


a) oj = — 1, x = —5C3X 6 — 45C2^ 4 — 15Cix 2 
+ ’ I ^12^32 < 6 

27 


b) ^ = 0, x 2 = 


V’32 > 0 V>32 < 0 

or < 27 

■ 012^32 <6 [ 6 < - 012^32 < — 


F_, ^32 < 0, ^12^32 < 


Li + L2 

Substitution Hi = x into equation 
of motion and constraint leads to 


a) cj = —1, 

W+, 


x 2 = 


Cl <0 

C2>0 


X < 


5C? (81) 

4C2 


, Cl < 0 or 

V>12 


Cl >0 

C 2 <0 


b) x = 0 , x 2 =-—, 4 >i 2 < 0 


Similar to the (7+l)-dimensional case, one can see that L 1 + L 2 +L 3 case is more abundant with solutions. 
We summarize all solutions obtained in Table IV and will discuss the differences and familiarities in detail 
in the appropriate section. 


Table IV. Summary of (6+l)-dimensional solutions in Li + L 2 + L 3 and EGB. 


Splitting 

nonvacuum L\ + +2 + 

vacuum +i + +2 + T3 nonvacuum EGB 

vacuum EGB 

(2 - 2 + 1 + 1) (68) 

yes 

no 

(2 + 2 + 2) (70) 

yes 

no 

(3 + 2 + 1) (72) 

yes 

no 

(4 + 1 + 1) (73) 

yes 

no 

(3 + 3) (75) 

yes 

(4 + 2) (77) 

yes 

(5 + 1) (79) 

yes 

no 

(6 + 0) (81) 

yes 


VI. SPACES WITH THREE-DIMENSIONAL ISOTROPIC SUBSPACE 

Now we want to turn our attention to the special case with three-dimensional isotropic subspace. Despite 
the fact that in the previous section we gave complete description of all possible cases, due to importance of 
the case with three-dimensional isotropic subspace - it could give a rise to successful dynamical compactifi- 
cation scenario where three dimensions expand isotropically while the remaining are contracting, and they 
can contract anisotropically - we decided to devote a separate section for this comparison. 

In (7+l)-dimensional L\ + L -2 + L 3 case we can have above mentioned scenario as a special case of (52) 
when one of the { y , z, u\ is equal to x > 0 (turning it into (3 — 2 + 1 + 1)) and the remaining two are 
negative. For (3 — 3 + 1) case (54) we require x > 0 with y < 0 and with no additional constraints retrieve 
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dynamical compactification of all but three space directions. Next case - (3 + 3 + 1) (55) - once again could 
be easily compactified - for nonvacuum solution we require x > 0 with z < 0 (we can always do this since 
z € R); we also need y < 0 which could be achieved with 032 > 0 (see (55)) and so 012032 < 1. For vacuum 
solution one can see that both branches for £ satisfy dynamical compactification as long as z < 0. 

The case (3 + 2 + 2) (57) is more interesting - indeed, for nonvacuum case we need B i(x) ^ 0 while 
y, z < 0 for x > 0; analytical solution of this inequality is next to impossible to obtain so we enumerate 
possibilities over a reasonable range of 0i2, 032 and x > 0. Our analysis suggests that there is a upper 
bound on x: x < 6.5 1 and allowed regions for 0i2 and '032 depends on x. One can simply verify that vacuum 
Li + L2 + T3 solution is also fine. Nonvacuum EGB solution is fine as long as y, z < 0 - constraints on both 
could be easily derived from (57). 

Finally, (4 + 3) case (60), and now we require x < 0 while y > 0. Similarly to the previous case we 
numerically investigated allowed regions on {032, 012, x} space, which lead to compactification. We found 
that for any x < 0 there exists region on {032, 012} plane which lead to y > 0, so that the limit similar to 
the previous case does not exist in this one. To support our claim we rewrite (60) for a and f3 (keep in mind 
(25)): 


a = — 


9 y 3 + 22>y 2 x + 27 yx 2 + 12x 3 


P = 


64xy 4 + 204x 2 y 3 + 264x 3 y 2 + 144x 4 y + 25x 5 ’ 
1 5x 4 + 12 x 3 y + llx 2 y 2 + 6 xy 3 + y 4 


(82) 


x, y € R. 


24 yx 3 ( 6x 4 + 36x 3 y + 66 x 2 y 3 + 51xy 3 + 16y 4 ) 

As for EGB case, one can easily see from (60) that both vacuum and A-term case always have successful 
dynamical compactification: for nonvacuum with x < 0 we can choose branch y = —3x + \/6x 2 — 0i2 to 
ensure y > 0 while for vacuum both £ give opposite signs to x and y. 

We summarize all possible solutions in Table V. One can clearly see that the situation with solutions 
which allow “good” compactification follows the general trend - we have L\ + L 2 + T3 case more abundant 
on the solutions; as of the comparison of A-term solutions with vacuum, the former are more abundant then 
the latter. 

Now let us perform the same analysis for (6+l)-dimensional solutions. The case with (2 — 2 + 1 + 1) 
splitting has only nonvacuum T1 + L2 + T3 case (68) and it could be brought to (3 — 2 + 1) with special choice 
of {y, z} = x > 0, while the remaining Hubble parameter is negative. Nonvacuum L\ + L 2 + T3 (3 + 2 + 1) 
case (72) requires additional constraint coming from y < 0 requirement plus we require z < 0; vacuum 
L1 + L2 + T3 and nonvacuum EGB cases satisfied by default: the former of them has the only solution which 
already has opposite signs for x and y (plus 2 < 0) while for the latter we just need z < 0. Final case to 


1 Keep (26) in mind for physical value of x. 
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Table V. Summary of (7+l)-dimensional solutions which could lead to successful dynamical compactification in 


L\ T L 2 + T3 and EGB. 


Splitting 

nonvacuunr Li + L 2 + L 3 

vacuum Li + L 2 + L 3 

nonvacuum EGB 

vacuum EGB 

(3-2 + 1 + 1) (52) 

special choice, say, 
y = x > 0, z, u < 0 

no 

(3-3 + 1) (54) 

y < 0 

no 

(3 + 3 + 1) (55) 

Z < 0, 032 > 0, 012032 < 1 

z <0 

no 

(3 + 2 + 2) (57) 

x < 6.5 (see test for details) 

yes 

constraints from y < 0 

no 

(4 + 3) (60) 

yes (see text for details) 


consider is (3 + 3) (75) and there are four subcases: nonvacuum L\ + L 2 + L 3 - using enumeration technique 
in reasonable range for x > 0, 0i2 and 032 we proved that for any x > 0 there exist area in {'0121 032} space 
which corresponds to y < 0; additionally we can, similarly to (4 + 3) case, rewrite (75) for a and j3 (again, 
keep in mind (25)): 


_ 3x 2 + 4 xy + 3 y 2 _ 1 x 4 + 4x 3 y + 5 x 2 y 2 + 4 xy 3 + y A , . 

a 16x 3 y + 28x 2 y 2 + 16xy 3 ’ 24 x 3 y 3 (4x 2 + 7xy + Ay 2 ) 

Vacuum L\ + L 2 + L3 case - from P- 2 (x , y) = 0 with use of y = £x substitution we get £ € [—1-5; 0) so 
that we also have compactification; nonvacuum EGB: y = — 2x — yj 3x 2 — 0i2 branch always gives 11s y < (1 
for x > 0; and finally vacuum EGB case also does possess successful compactification since y = £x and 
£ = —(3+ \/5/2 < 0) both. So that all (3 + 3) cases have dynamical compactification. Again we summarize 
all possible solutions in Table VI. 


Table VI. Summary of (6+l)-dinrensional solutions which could lead to successful dynamical compactification in 

L\ + 1*2 + T3 and EGB. 


Splitting 

nonvacuum L\ + L 2 + L 3 

vacuum L\ + L2 + L 3 nonvacuum EGB 

vacuum EGB 

(3-2 + 1) (68) 

special choice, say, y = x > 0 , z < 0 

no 

(3 + 2 + 1) (72) 

constraints from y < 0 + z < 0 

z <0 

no 

(3 + 3) (75) 

yes (see text for details) 


VII. CONCLUSIONS 

In this paper we studied higher-dimensional exponential solutions for EGB and Lovelock gravity with 
third correction. As we mentioned in the Introduction, both cases - EGB and mode general Lovelock gravity 
- are eligible and well-motivated, and they demonstrate different properties of the solutions. We developed 
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a scheme which allows one to find all possible exponential solutions in any number of dimensions and with 
all possible Lovelock terms. The scheme was tested on (6+1)- and (7+l)-dimensional space-times and all 
possible solutions for them were found - they are summarized in Tables III and IV. We can see that these 
solutions are more abundant in the presence of all three Lovelock terms as well as in the presence of non¬ 
zero cosmological constant in the action. In the most general case in (7+l)-dimensions we have 10 possible 
splitting of initial 7-dim space into different isotropic subspaces (which number does not exceed 5), a vacuum 
general third order Lovelock theory has 7 possibilities (with no more than 3 isotropic subspaces), non-vacuum 
EGB theory has 6 cases, and a vacuum EGB - only 3 of them (the latter theory has solutions only with 
two different subspaces). The same tendency appears in the ( 6 +l)-dimensional case with 8 solutions in the 
general non-vacuum case and 3 solutions in the vacuum EGB theory. 

As a special case we separate the solutions which allow three-dimensional isotropic subspace with different 
signs for Hubble parameters corresponding to this subspace and all remaining directions - this requirement 
ensures that three-dimensional subspace will expand (given positive sign for the corresponding Hubble 
parameter) while in the remaining directions Universe will contract. In turn, this allows us to consider this 
situation as successful example of dynamical compactification - indeed, we have three dimensions expanding 
and the remaining contracting - though, they could contract anisotropically. We summarized our findings 
of the solutions of this type in Tables V and VI. Similar to general solutions, the more terms in the action 
of the theory we consider, the more possibilities for compactification we have. 

This finalize first part of our study of exponential solutions in Lovelock gravity. By now we completed 
investigation of exponential solutions in EGB gravity, described both non- and constant-volume solutions. 
In this paper we generalized non-constant-volume solutions on the case of Lovelock gravity. By comparing 
solutions and their abundance for EGB and L\ + L 2 + L 3 one clearly see that more general Lovelock case 
has more different solutions, which is partially due to increased number of degrees of freedom (through 
increased number of couplings), and with growth of the number of dimensions the highest possible Lovelock 
order is also increasing, resulting in growing number of possible solutions. This makes the same solution 
with three-dimensional isotropic subspace more “probable” to appear in case with more Lovelock corrections 
taken into account, which, in turn, could be used as kind of anthropic argumentation for consideration of 
Lovelock and other higher-order gravity models. 
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